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We analyze the generic structure of Einstein tensor projected onto a 2-D spacelike surface S 
denned by a unit timelike and spacelike vectors u and n respectively, which describe an accelerated 
observer (see text). Assuming that flow along u defines an approximate Killing vector £, we then 
show that near the corresponding Rindler horizon, the flux j a — G^£, b along the ingoing null geodesies 
k, i.e., j • k, has a natural thermodynamic interpretation. Moreover, change in cross-sectional area 
of the k congruence yields the required change in area of 5 under virtual displacements normal to it. 
The main aim of this note is to clearly demonstrate how, and why, the content of Einstein equations 
under such horizon deformations, originally pointed out by Padmanabhan, is essentially different 
from the result of Jacobson, who employed the so called Clausius relation in an attempt to derive 
Einstein equations from such a Clausius relation. More specifically, we show how a very specific 
geometric term [reminiscent of Hawking's quasi-local expression for energy of spheres] corresponding 

to change in gravitational energy arises inevitably in the first law: diJc/dA oc f d 2 x^/a (see 

H 

text) - the contribution of this purely geometric term would be missed in attempts to obtain area 
(and hence entropy) change by integrating the Raychaudhuri equation. 

PACS numbers: 04.62.+v,04.60.-m, 04.70. Dy 



I. INTRODUCTION 

The study of thermodynamic aspects of black holes over the past decades has given several insights into the 
nature of gravity as described by Einstein's General Relativity, and is expected to be a crucial link in constructing 
a quantum theory of gravity (see [1] for a recent review and references). In a paper more than a decade back [2], 
Jacobson speculated that it might be possible to invert the logic of the "physical process" version of the laws of 
black hole mechanics, developed by Wald, and by applying it to local Rindler horizons, one can derive Einstein field 
equations from Clausius relation, TdS — &Em-, where Em is related to matter flux (and vanishes when T a t = 0). 
The essential new idea introduced by Jacobson was that of local Rindler horizons in a small patch of spacetime which 
can be approximated as flat once one has set the acceleration length scale appropriately. (See Appendix A 1 for an 
elaboration on this construction.) Einstein equations would then emerge as consistency conditions on the background. 

In a later paper [3] , Padmanabhan pointed out that if one actually looks at the structure of Einstein tensor near 
a spherically symmetric horizon, it has the form TdS" = dEc + PdV, where Eq is associated with horizon energy 
(and unlike Em, Eq ^ when T a b = 0) and P with matter flux (these are defined below). In fact, the above relation 
has been shown to hold for a wide class of horizons, including arbitrary static horizons in Lanczos-Lovelock theory 
as well. This result looks different from what Jacobson had started with to deduce the null-null part of Einstein 
equations - specifically, the energy term Eq has nothing to do with Em, which is more like the PdV term but with 
different interpretation in terms of matter flux. So, while the Clausius relation seems to yield null-null component of 
Einstein equations, the Einstein tensor itself has a very different structure. It is important to relate these results and 
understand where the difference comes from, which we intend to do in this note. 

Before proceeding, we would like to clarify an important point so as to put the analysis presented here in proper 
perspective. To begin with, we must mention that our main emphasis here is not to analyze pros and cons of one 
method over the other, but rather to clarify why they differ and to characterise the diffcrence(s) from a physical point of 
view. It is indeed true that a priori there is a difference between the approaches of Jacobson and Padmanabhan; while 
Jacobson's analysis concerns deriving Einstein equations from Clausius relation, Padmanabhan's result demonstrates 
that Einstein equations on horizon is same as the first law of thermodynamics. However, once the physical content 
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of Einstein equations has been claimed to be equivalent to a particular thermodynamic relation, one would have 
expected a mapping between the two results, unless there are subtle differences at a fundamental level. Indeed, the 
TdS term in the thermodynamic relation is fairly unambiguous, so that the remaining terms in the equations must 
correspond in some manner. If they do not, then it implies that there is a difference at a conceptual level, which is 
what we shall show in this note. We shall show that the difference arises in the particular manner in which matter 
fluxes across the horizon are treated. Specifically, Padmanabhan's result arises due to deformations of the future 
horizon normal to itself, generated by ingoing null geodesic congruences, and this yields the force term PdV in the 
final result. We discuss in some detail the resulting difference in physical interpretations. Furthermore, we show that 
the additional term dEc is essentially the change in quasi-local energy associated with the horizon 2-surfacc, and is 
related to horizon topology; more precisely, we show that dE G /dX oc J d 2 x y/a ^R. To summarize, we shall do the 
following in this note: 

1. Clarify the role of horizon deformations to be considered in a Rindler patch when matter crosses the future 
Rindlcr horizon of the observer. 

2. Clarify the differences between the "heat flux" term of Jacobson, and the "PdV"' term of Padmanabhan, and 
highlight the physical implications. 

3. Indicate clearly that the change in area of a horizon cross-section is determined by the expansion (and not its 
first derivative) of the ingoing null congruence normalised to have unit Killing energy. 

4. Give an explicit expression for the expansion 9 for the congruence mentioned in the previous point, in terms of 
combination of curvature tensor components (see Eq. (5) below), and compare with the corresponding combi- 
nation occuring in the Raychaudhuri equation. In particular, the area change involves not just the Ricci tensor, 
but also the Riemann tensor - a point which is of relevance in the context of deriving field equations from 
thermodynamics. 

5. Exhibit explicitly the "thermodynamic" structure of Einstein tensor and show that there is a term corresponding 
to quasi-local energy of the horizon, which must be separately accounted for when considering energy flow across 
the horizon. 

6. Show that, when using the Raychaudhuri equation with our prescribed null congruence, the O(A) term does not 
vindicate or necessitate setting the expansion to zero. 

We shall address all the above points in the following sections. To avoid distraction from the main points, we have 
relegated most of the mathematical details to appendices. Before proceeding, let us also clarify the restrictions on 
the local frame of the accelerated observer that we shall impose. The most important restriction is that of staticity; 
that is, we shall assume that, in the local coordinates near the observer worldline, one can define an approximate 
timelike Killing vector field. Consequently, the near horizon geometry is assumed to be static. For static spacetimes, 
we shall use, for the near horizon metric, the form: ds 2 = — N 2 dt 2 + dz 2 + aABdy A dy B , with the Taylor expansions 
for N and gab derived by Visscr ct al [4]. As our discussion will make clear, the above form of metric is just a good, 
convenient parametrization - the final results are of course stated in a manifestly tensorial form. The only crucial 
input is staticity, which requires a satisfactory notion of a timelikc Killing vector which is hypersurfacc orthogonal, 
and a spacelike surface whose unit normal points in the direction of acceleration. 

II. THE NULL BASIS NEAR A HORIZON 

Let us concentrate on the future horizon H. of the right Rindlcr wedge, which is generated by outgoing null rays. 
The most natural transverse null vector for H is therefore defined by affinely parametrized [5] ingoing null geodesies, 
fc, and can be chosen to be: k = N^ 1 ^ — n). Here, N = \J — £ 2 , u = £/N, and existence of a local timelike Killing 
field £ (which generates local Lorentz boosts) is assumed. Also, n is the unit normal in the direction of acceleration of 
u. The existence of £ is also assumed in the work of Jacobson [2], and without this no further progress can be made. 

The choice of normalization is such that k ■ £ = —1, implying that k has unit Killing energy. It must be also 
noted that the corresponding outgoing null rays are given by I = A^ _1 (tt + n); we note that, N 2 l — > £ on %. More 
precisely, N 2 l become tangent to horizon generators [the vector I itself does not, since I •£ = — 1 by construction]. The 
standard Rindlcr transformations in the local inertial fram (LIF) has an additional parameter re which characterizes 
orbits of Lorentz boosts and generates constant acceleration trajectories. In inertial coordinates (T, X, Y, Z), we have 
k = k (X + T) -1 (dr — dx) —> (2reX) _1 (dr — dx) on W, i.e., T = X. In fact, we could as well have used I below 
for discussion without any change in the final result, but this would be a weird thing to do since these geodesies 
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behave badly near the future horizon. Specifically, I = k~ 1 {X — T)^ 1 (8t + dx), so the components in the locally 
inertial coordinates blow up at X = T. 

We shall now demonstrate that the expansion of fc (or i) governs the changes in cross sectional area of the horizon. 
To exhibit the result for both fc and I simultaneously, we write fc e = N^ 1 (it — en), where e = +1 corresponds to fc 
and e = —1 to I. First, note that, in terms of the covariant derivative ^ D compatible with t ^constant hypersurface, 
we have 

{ ^D-n = {g ab + u a u b )^ a n b 

= V • n — n ■ a (1) 

Now evaluate 

V-fc £ = -eJV _1 V -n + Nk.-VN^ 1 

1 (o\ an k e ■ V7V 

= — e — y, D-n—E 

N N N 

= - e — &D-n (2) 

N V ; 

where we have noted that a = VN/N, so that the last two terms in second equality cancel. Now, since the t =constant 
metric is dz 2 + a ab&V A dy B ', the n = d z congruence is an affinely parametrized geodesic congruence, and therefore 
we can use the standard interpretation of ^'D ■ n in terms of fractional rate of change in "volume" of z =constant 
surfaces, which corresponds to the 2-D manifold described by the metric gab- Hence, we finally get 

V • fc c = e-—d z In \fo = —€—— In \fa (3) 
TV dA 

which is the desired result. (We have used Ndz = dA in arriving at the second equality, see Appendix A 2 for details.) 
This straightforward evaluation should leave no doubt as to how the change in cross sectional area of the horizon is 
actually described by considering ingoing (or outgoing) null geodesic congruences, constructed in the manner we have 
described. In fact, the choice of ingoing congruence fc is also strengthened by some old results due to T. Dray and G. 
't Hooft [6] , which clearly shows that a massless particle falling into a Schwarzschild black hole corresponds to a shift 
in the ingoing Kruskal coordinate, the shift being proportional to the particle energy. 



III. THE THERMODYNAMIC STRUCTURE OF EINSTEIN TENSOR 



We shall now analyze the near-horizon form of Einstein tensor, and reveal how its thermodynamic structure emerges. 
Before proceeding, however, we wish to clarify an important point concerning the variations we shall be considering. 
We shall base our discussion on the ingoing null geodesies fc of the previous section, satisfying fc • £ = — f . As should be 
evident from comments in the previous section, the entire analysis can be repeated in a straightforward manner using 
outgoing null geodesies I satisfying I ■ £ = — 1; the only difference is the change in sign in n at various intermediate 
steps, while the final result remains unchanged. The reason for using ingoing null geodesies fc, as mentioned above, is 
that these have components which are well behaved at the future horizon in the locally inertial coordinates; the only 
crucial thing is the normalization based on unit Killing energy. 

We begin with the following (exact) identity [see Appendix B for a proof]: 

G ab g ±ab = 2 (R ab Ck b - R abcd u a n b u c n d ) - 

- N 2 R ab k a k b + U[K, fc] (4) 

where g^ b = —u a ii b + n a n b is the metric on the surface orthogonal to the horizon, and H[K, fc] = /(fc) — f(K) — 
<j>(K), with f(K) = K 2 - K 2 U (similarly for /(fc)), and <j>{K) = n^ n P (K»K vp - KK W ). Here, and k AB are 
extrinsic curvatures of level surfaces of u embedded in 4-D spacetime, and of n embedded in the resultant 3-D space, 
respectively Note that the above expression is true for an arbitrary spacetime without any geometric constraints 
imposed so far. 1 



1 In particular, for a flat 3-D space in a flat 4-D spacetime, one obtains ( 2 ' R = f(k), which is essentially the content of Gauss's Theorema 
Egregium. 
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We shall now impose the condition of staticity, that is, we shall require that the near horizon geometry, to a sufficient 
approximation, has a local timelike Killing vector field. In that case, we can show that (see Appendix C), on the 
horizon z — > 0: 

RabCk b - R ab cdu a n b u c n d = k^- In yfd (5) 

dA 

The above equation gives the derivative of area (rather than it's second derivative) in terms of curvature components, 
and deserves several comments, which we list below: 

• It clearly shows that the change in cross sectional area (obtained by integrating \fa over transverse coordinates) 
of the k (or the I) congruence (normalized so as to have unit Killing energy), on a cross section of H, depends 
on a very different combination of Riemann tensor components than the one occurring in the Raychaudhuri 
equation [which only involves Ricci tensor, R a i,k a k b ]. 

• Raychaudhuri equation gives second derivative of area and our analysis above shows that "integrating" it naively 
to obtain the first derivative will, in general, be tricky. Indeed, the null-null component does not appear in the 
above equation at all! In section IV, we shall present an analysis a la Jacobson using Raychaudhuri equation, 
which should clarify further what is going on here. 

(This and the previous comment arc important particularly when we consider Jacobson's argument and compare 
it with our result, see section IV.) 

• The appearance of R a bcdU a n b u c n d also must be highlighted; one could have simply ignored this term by de- 
manding it to be small, and calling this demand a further restriction on the definition of a local Rindler horizon. 
This, however, would be adhoc, since for Schwarzschild horizon, it involves d 2 (l — 2M/r). Indeed, as is evident 
from above, there is actually no need to throw away this term, since it occurs in just the right combination in 
Einstein tensor so as to give the change in area correctly. 

• Even if we did throw away the R a bcdU a n b u c n d term, we are left with R a b£ a k b which has nothing to do with the 
null- null component of Ricci [recall that fe ■ £ = — 1] . 

Proceeding to the main analyis, note that if R a bk a k b [and hence G a bk a k b ] is finite on the horizon, then the cor- 
responding term on RHS of Eq. (4) is 0(z 2 ). Also, H[K,k] is ignorable because it is 0{z 2 ). This comes about as 
follows: K^ v is zero due to staticity. On the other hand, kAB cx d z a ab is O(z) since, from the Taylor expansion of 
area, gab = (^-independent part) +0(z 2 ) (see Ref. [4]). Since H[K, k] is quadratic in Aus, it is 0(z 2 ). So we finally 
obtain: 



where we have defined P = (l/2)T ab g ±ab . The differential version of the above equation (multiplying it by dA) yields 
Padmanabhan's result: 

PdV = TdS - dE G (7) 

Having established the above relation, we can ask how general it is. It might seem that the result is very specific to 
Einstein gravity since in arriving at it, we used Eq. (5) for change of area, and in Einstein gravity horizon entropy is 
proportional to area. We could therefore relate entropy change to area change and derive the result. However, when 
one goes beyond Einstein theory, entropy is no longer proportional to area but is instead given by Wald entropy. It 
is therefore quite a non-trivial fact that exactly the same result can be proved for a much larger class of lagrangians 
- the so called Lanczos-Lovelock (LL) lagrangians - for which the horizon entropy is given by a non-trivial function 
of area. Once again, we find that the near-horizon structure of field equations for LL actions can be cast in the form: 

PdV = TdS LL - dE (G)LL (8) 

and the resulting expressions for Sll and -E(g)ll turn out to be [7]: 

Sll oc j ^4?-i 2) (9) 
(d£ G /dA) LL cx j V^L^-V (10) 
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We see that S is precisely the Wald entropy, whereas Eq gives the correct expression for quasilocal energy when 
applied to known black hole solutions. 2 

Before turning to Raychaudhuri equation, let us make another relevant comment: It is easy to sec, from our 
definitions, that G a b9 ±ab = —2G a b£, a k b + N 2 G a bk a k b . Therefore, provided G a bk a k b is finite on the horizon, one 
obtains, in the limit TV — > 0: G a b^ a k b — > — (l/2)G a bg ±ab (which, incidentally, is the so called work function, W, 
defined by Hayward in the context of spherically symmetric, dynamical horizons [8]; note that we have not assumed 
spherical symmetry to obtain Eq. (6)). On the horizon, we therefore have natural interpretation for this term as the 
force acting on the horizon in the direction defined by fc. Let us also mention its form for an ideal fluid, described by 
Tab = PoVaVb + Poidab + v a Vb), where v a is the fluid 4- velocity, and we assume for simplicity that it lies only in the 
u-n plain. Then, a trivial calculation shows that T ab u a u b = j^eiiPo + Po v rei) an d T ab n a n b = 7 2 e ;(Po + Po v rei)i wnere 
Irei = —u ■ v - We then immediately obtain P = (l/2)T ab g ±ab = (p a — po)/2. 

It is also instructive to compare this analysis with the one given by Jacobson, in which case the most natural 
starting point would be the Raychaudhuri equation. We do this in section IV. We shall show that, for our k (or I) 
congruence, the starting assumption of equating TdS with matter flux gives, at O(A ), a relation which is inconsistent 
with the algebraic identity obtained in this section. However, if one makes further approximations and ignore certain 
terms, then we do recover the null- null part of Einstein equations at O(A), although in a manner completely different 
from Jacobson's, since our analysis is not based on the null generators. Most importantly, we do not require the 
vanishing of expansion of the null congruence at all. Before proceeding, we must emphasize that, in the next section, 
we shall be trying to follow Jacobson's reasoning in our setup; the final results and implications must, of course, be 
interpreted keeping this in mind. Needless to say, our main emphasis is towards trying to understand why there are 
differences between the work and energy terms in the two approaches; the answer, as we hope this note would make 
evident, lies in different ways of treating fluxes across the horizon. 



IV. ANALYSIS BASED ON RAYCHAUDHURI EQUATION 



In this section, we turn to Raychaudhuri equation, in an attempt to understand better the difference between above 
result and Jacobson's derivation of the null-null component of the field equations. To do so, we repeat Jacobson's 
analysis using the k congruence; this should indicate where the difference lies. Once again, it is worth emphasizing 
that we will obtain the same results upon using the outgoing I congruence of unit Killing energy. As we have shown 
above, the Einstein tensor on the whole has a much richer structure due to the presence of the ^ R term, which we 
would want to explore further. Unfortunately, Raychaudhuri equation, as we will see, has nothing much to say about 
this term, but our analysis will shed some light on the role of certain assumptions in Jacobson's derivation, and also 
the differences between the work term as well as horizon energy. 

Start with the equation defining variation of area in terms of expansion 9 of a congruence of ingoing null geodesies 
w.r.t. the affinc parameter A along k (see Appendix A 2 for more details). Assuming that entropy is proportional to 
area, this gives: 

T H dS = a' 1 J 6»d£ dA (11) 

where a = (8ttcL 2 /H)/k, and the integration is over the null 3-surface generated by the cross-section of a bundle of 
ingoing null geodesies k across an afline distance A. The horizon is at A = 0. Now expand 9 

9(X) = 6(0) + 9{0)X + i(9(0)A 2 + 0(A 3 ) (12) 

in obvious notation. Now we can use Raychaudhri equation to substitute for the first derivative of 9 evaluated at 
A = 0. That is, 

9(0) = -^(0)-[R ab k a k b ] x=0 (13) 
where we have ignored shear and rotation for the time being (which is also an assumption in Jacobson's work). 



2 A general definition for quasilocal energy, such as Hawking's definition for Einstein theory, is not available for the LL actions; in fact, 
ours can be taken as a natural generalization of Hawking's quasilocal energy for LL actions. 
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Now consider the heat flux through d£ dA: 



dQ 



T ab Ck b dS dA 



for which a similar expansion gives: 



dA 



0(A 2 ) 



(14) 



(15) 



A=0 



Following Jacobson, we now impose the Clausius relation, TdS = dQ and equate equal powers of A on both sides. 
That is, 



0(O)+0(O)A 
This gives 



-d{0)X 2 + 0(A 3 ) 



dS dA 



[T ab Ck b ] x=Q + X 



T ab Ck b 



dA 



0(A 2 



A=0 



0(A°) 
0(A X ) 



e(o) = a[T ab Ck b ] x=0 



--9 2 (0)-[R ab k a k b ] x= =a 



T ab Ck b 



dA 



dSdA (16) 

(17) 
(18) 



A=0 



Using Eq. (17) to replace 2 (0), this becomes 



o(a x ; 



-- [aT ab ek b ] 2 x=0 - [Rabk a k b ] x=Q = a 



jrT ab ek b 
dA 



A=0 



The relevant points to note here are: 



- On the horizon, nXk a goes to £ a , that is [«Afc a ] A=0 = £ a (see Appendix A 2) , which is obviously a O(X ) 
expression and NOT O(X). This is a key difference from Jacobson's argument, arising because Jacobson considers 
fluxes along generators k a of the horizon. [In that case, nXk a ~ £ a is valid all across the Killing horizon (A being 
the affine parameter along the generators k a ). This then necessitates that expansion of the generators vanish at 
the bifurcation surface A ~ (corresponding to T = = A), since the matter flux term becomes O(A).] In our 
opinion, since one would expect to associate entropy with cross sections of arbitrary null vectors in an arbitrary 
curved spacetime, such an assumption on the expansion is restrictive. 

In our setup, it would actually be incorrect to deduce that [T a b£, a k b ] X Q is O(A); as seen from above, this term 
is in fact related to 0(0), which in general does not vanish. In fact, one would expect arbitrary null congruences 
to block information of a certain region of spacetime from a class of observers; for such congruences, there is 
actually no need to constrain the expansion to vanish. 

So, whether Einstein equations come out at O(A) depends on [d\ (T a b^ a k b )] X _ Q = [A; a V a (T a b£ a £; b )] x . In general, 
it is not at all obvious what this will lead to, but let us consider this term in more detail: 



fc c V c [T ab £ a k b 



ek b k c V^ab + T ab Ck c V c k b + T ab k b k c V c C 

ignore —0 



(19) 



where we have ignored the derivatives of T ab , which is justified in the approximation in which we are working here. For 
consistency, one must then also ignore the T 2 b term on the LHS of Eq. (18). Now concentrate on the term involving 
k c V c £, a , which is to be evaluated at A = after computing the derivative. This term at A = can be shown to give 
— K,k a . Therefore we obtain: 



fc c V c [T ab Ck b ] « -KT ab k a k b 



(20) 



The last term in [d (T ab ^ a k b ) /dA] therefore reproduces precisely the contribution which would come from calling 
(in our case incorrectly) [T a b^ a k b ] x _ as a 0(A) term! Therefore the 0(X 1 ) becomes [Note that a = (8ncL 2 / K) / k]: 
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OiX 1 ) : -[R ab k a k b ] x=0 = -(87rcL 2 p /h)[T ab k a k b ] x=0 

thereby giving the null-null component of Einstein equations yet again, as Jacobson had obtained, but in a completely 
different manner! Also, note that the above relation is applicable all along the future horizon, and not just near the 
bifurcation point, in so far as the notion of local, static Rindler horizon remains well defined. 
The discussion above clearly implies that: 

1. In our setup, Einstein equations do NOT necessarily follow from TdS — d-E^atter; this is so since one could also 
have included an additional term, involving curvature tensor, which will only modify the O(X ) term, which is 
anyway ignored while evaluating the Oi^X 1 ) contribution. In fact, as we demonstrated in the previous section, 
such a term is present in Einstein equations, and corresponds to a change in "gravitational energy" . 

2. The above claims can be explicitly demonstrated by evaluating Einstein equations near a static horizon, in 
which case a very natural energy term is picked out; the resultant equation, in fact, can be thought of as 
TdS* — diShorizon = P&V = -dBmattef Note that the O(X ) piece above can be rewritten as (fi;/27r)d(<L4./4) = 
— T a b£, a k b 5A, which is, of course, Padmanabhan's result without the (*)R term. Of course, since the latter result 
is algebraically correct (as we showed in the previous section), this actually makes the O(X ) contribution (and 
by implication, the starting relation) incorrect at a purely algebraic level - one needs the ^ R contribution for 
consistency [the minus sign above is due to k being an ingoing congruence; as an interesting aside, let us also 
mention that this fact is one of the "boundary conditions" for dynamical/isolated horizons imposed by Ashtekar 
et. al.]. 

In fact, we feel that the " ^R" term reinforces the well known quasi- local character of gravitational energy; an 
ultra-local description of energy balance might therefore be a bit tricky. 

3. Let us also highlight the role of the ingoing congruence k. The other null congruence N 2 l generates the horizon, 
and since we have taken great pains to construct a local Killing horizon, the expansion along the generators would 
vanish everywhere on the horizon so long as our local constructs make sense; when they don't, we cannot even 
talk about local Rindler observers and local Killing horizons. Instead, the congruence k captures information 
infinitcsimally away from the horizon in the direction normal to it [that it is a natural "normal" is easily seen in 
inertial coordinates, where future horizon is T — X — whose normal is clearly along fe.], and provides a natural 
flow to define variations of various quantities. We hope to have shown that it is this congruence which gives, in 
a sensible manner, the change in a cross-section of the horizon, as well as the matter flux across it. Of course, 
the result can also be justified by applying to event horizons of a black hole solution of Einstein field equations. 
For a stationary black hole horizon, the only sensible change in area when matter crosses the horizon is given 
by expansion of fc, since the expansion of horizon generators vanish for a Killing horizon. Moreover, as we have 
already mentioned before, it has been shown rigorously by T. Dray and G. 't Hooft that a massless particle 
falling into a Schwarzschild black hole corresponds to a shift in the ingoing Kruskal coordinate, the shift being 
proportional to the particle energy. This further strenghtens our motivation for using the ingoing congruence 
for horizon deformations. 

4. One of our main conclusions, as far as comparison with Jacobson's analysis is concerned, is the following: The 
difference between Einstein equations being identical to the first law of thermodynamics, as was pointed out by 
Padmanabhan, differs from the Clausius relation of Jacobson due to difference in the manner matter fluxes 
across the horizon are defined. Our analysis seems to be closer to the one in Refs. [8] and [12]. 

5. There are also other significant issues which go beyond the algebraic ones we have mostly concentrated on till 
now. In using the Clausius relation, one is trying to derive field equations from a starting thermodynamic 
relation. In such a case, one has to change the starting relation depending on what type of congruence one 
chooses. These additional terms are then interpreted as dissipation terms and are accounted for by adding 
suitable entropy production terms in the Clausius relation. However, our analysis above shows that, for a 
suitably defined horizon (with static near-horizon geometry) , the field equations take the form of the first law 
of thermodynamics without involving any additional terms, under prescribed horizon deformations. The only 
sensible quantity to concentrate on, while comparing these two approaches, is the TdS term (which can have no 
ambiguity once the entropy density is suitably defined, and which can be verified by applying it to known cases 
of black hole horizons); in our case, this term is derived using the expansion 9 alone, rather than its derivative, 
as is required while using Clausius relation. However, as has already been emphasized several times before, it 
must be remembered that the difference originates in using different definitions for matter fluxes (and not due to 
any in-correctness in any of the two approaces); the physical motivations for the choice used here are mentioned 
in point (3) above. 
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Before moving on, we would like to point out that there have been some recent attempts to derive field equations 
for higher derivative gravity theories [13] along the lines of Jacobson. In these works (except Padmanabhan's, see 
below), the starting point is Clausius relation along with Wald's definition of entropy in terms of Noether charge of 
diffcomorphism invariance. However, while the definition of matter flux is similar to Jacobson's, the Raychaudhuri 
equation is never invoked, thereby avoiding any need for assumptions such as vanishing expansion etc. present in 
original Jacobson work. It would be worth investigating further how the comments in this note are to be considered 
in the context of these recent attempts. (For one thing, the difference in the definition of matter flux remains.) We 
must, however, point out that the situation is far from clear since these works do not all agree with each other. 
For e.g., Parikh and Sarkar have pointed out issues with the Brustein and Hadad paper, whereas Padmanabhan has 
highlighted certain conceptual issues regarding interpretation of both these papers. Specifically, Padmanabhan has 
stressed the subtleties in interpreting these results as derivation of field equations from thermodynamics; the result, he 
argues, is better viewed as interpretation of field equations as a thermodynamic relation. More importantly, Brustein 
and Hadad derive their result using a different form for Noether potential (used to define entropy) than the other two 
papers; this calls for a more detailed and critical look at the analysis therein. Moreover, the fact that they still derive 
the same result implies non- uniqueness of the analysis (which a quick look at these papers will confirm), whereas in 
Jacobson's case, once the assumptions are stated, the analysis is unique. Hence, the status of these results in the light 
of the original Jacobson calculation remains unclear; the derivations are not only very different from Jacobson's, they 
are not even similar to each other! Perhaps the most important point indicating why these analyses are conceptually 
different from Jacobson's is f(R) theory: whereas the above papers derive the field equations from a Clausius relation, 
Jacobson and collaborators needed to add extra terms to the Clausius relation in their earlier paper [14], to proceed 
with the analysis. We hope further work will clarify these issues. 



V. IMPLICATIONS 



In this brief paragraph, we would like to emphasize the need for the analysis done in this note. Whether Einstein 
equations are just equations of thermodynamics in disguise is a well motivated question, and wc do agree that the 
answer to this might be yes. The important point realized by Jacobson in [2] while addressing this question, was 
to introduce local Rindler frames in an arbitrary curved spacetime, and use the thermal aspects of corresponding 
horizons as probes of the background curvature. However, one needs to impose certain restrictions to proceed from 
there, and to put the result in a physically relevant context. The necessity of highlighting such restrictions goes hand 
in hand with identifying specific geometric quantities with (variation of) the thermodynamic variables. 

In this sense, as we have shown, the expression for change in entropy of a cross section of a static horizon is related to 
very specific components of the Riemann tensor (and is readily verified for known black hole solutions). Once wc agree 
on these algebraic identifications, Einstein equations take the nice form of the first law of thermodynamics, provided 
one attributes to a 2-surface an energy proportional to its intrinsic curvature. Therefore, Einstein equations resemble 
the first law of thermodynamics in this very specific form. Of course, we have also shown that using Raychaudhuri 
equation also yields, at a higher order (and after justifying ignoring certain terms), the null- null component of Einstein 
equations; the crucial point, however, is that the null-null component of Ricci itself has no clear meaning in terms of 
change of entropy. Looked at in this way, the null-null part of Einstein equations seem to be a secondary consequence 
of the first law itself. Of course, one can reinterpret Einstein equations as representing some sort of a Clausius relation; 
if one insists upon doing so, one must redefine the matter flux suitably. It is easy to see that such a definition of flux 
will involve the trace of the matter stress tensor, and will not be equivalent to the heat flux defined by Jacobson; see, 
for example, reference [8] which gives one such definition for dynamical horizons, but assuming spherically symmetry. 
In fact, this can be easily demonstrated using our Eq. (5): 

(5)MH-S^(*** + 5«0 

« ^ (r ab - \Tg a A k a e (21) 

where R± is defined in Appendix C. Wc have used the field equations in the second equality, and the approximation 
is obtained after ignoring the R± term. Moreover, inverting the logic and deriving field equations in this latter case 
is again not sraightforward. Wc hope to have clarified all the above issues in the present note. 

At a more conceptual level, one of the possible implications of this note is that it might be necessary to adopt a new 
starting point if one wants to establish Einstein theory in terms of thermodynamics, in which case the thermodynamic 
structure of the Einstein tensor would serve as the most important supporting evidence. This point of view, of course, 
also applies to the wider class of Lanczos-Lovelock lagrangians for which similar results hold. A survey of some of the 
attempts that have been made along these lines can be found in [1, 9]. 
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Appendix A: Some mathematical details 
1. Approximate local Killing vector generating Lorentz boosts 

The existence of an approximate boost generating Killing held in a locally flat patch of spacetime is a crucial 
assumption which goes into applying notions of thermodynamics to local acceleration horizons, and it is therefore 
important to make precise the details of such a construction, which we do below. Around an arbitary event V in the 
spacetime, one can construct Riemann normal coordinates y k (with y k (V) = 0) in terms of which the metric takes 
the form: g ab = jj j — (l/3)R a cbdy c y d + 0(y 3 dR). One now wishes to find an approximate Killing vector field in this 
patch of spacetime, which would generate approximate Killing boosts along some direction, say y . That is, we want 
an analog of Minkowski vector held £ M = ^[y 1 , y°, 0, 0] where k is acceleration of the chosen Rindler observer. This 
can be done by writing £ fe as a series in y k (starting, of course, at linear order) and then checking how much freedom 
one has in the choice of coefficients so as to satisfy Killing equations, S a b — V a £b + Vb^ a = 0. The analysis shows that 
£ k can be made Killing only upto cubic order, that is, S a b = 0(y 2 ). Incidentally, this also implies, upon using the 
identity (valid for arbitrary vector field): V c V Q £ fc = R ba cnC - (l/2)(V a S' fcc + V c S ab - V&S ca ), that Vt,V a £ c - R ca bn£, n 
is 0(y). Wc also have: = ~RU n - V a S h a , where 5 fc a = S£ - (l/2)S? ff £. 

2. Parametrization of null geodesies 

[The discussion below mostly follows Wald's book [11], section 6.4: The Kruskal Extension.] 

Consider the simpler form of the Rindler metric: ds 2 = — x 2 dt 2 + dx 2 . The null geodesies satisfy the equation: 
i = ±(loga;)', where dot denotes derivative w.r.t. affinc parameter. Therefore, t = ± logx+const. for the null 
geodesies. Further, due to staticity, — ko = x 2 t is constant along a geodesic, so that dA = Cx 2 dt = ztCxdx. Hence, 
A = ±C(x 2 /2) + D is the affinc parameter. This is the choice we have made. It is further easy to see that the null 
geodesies are given by k l = C~ 1 [l/x 2 ,±l/x). Further, note that C _1 is the Killing energy associated with the null 
geodesic; imposing k ■ £ = —1 fixes C = 1 - that is, the multiplicative factor in the choice of affine parameter is fixed 
by requiring the null geodesies to have unit Killing energy. 

Note that the above parametrization is meaningful only for ingoing geodesies near future horizon (or outgoing 
geodesies near the past horizon). An alternate parametrization is, of course, in terms of the Killing time t. It is easy 
to see that this leads to A cx exp (±2i)+const, which is perhaps a more familiar parametrization. Finally, it is easy to 
deduce the following limits for the contravariant and covariant components of the vector k: (a) [a; 2 ^ 1 ]^^ = [1,0] = 
(b) [xki} x=0 = [0, ±1] = rii. These limits clearly confirm that x = is a null surface, and also illustrates the subtlety 
associated with taking such limits in a coordinate system which is singular (that is, when the metric or its inverse 
blows up in some region, here x = 0). 

Appendix B: Gauss-Codazzi decomposition relations 

[For greater clarity and notational convenience, we use boldface subscripts for contraction on corresponding vectors 
in this and the following appendices; for e.g., R^k = Rab£ a k b .] 

Begin with the Gauss-Codazzi expression for the 4-D Ricciscalar R with respect to a foliation defined by a timelike 
unit vector u: 

R = 2e u R ab u a u b + ^R + e u [K^K^ - K 2 ) 

[See, for e.g., [10] sec 3.5.3, page 78.] Similarly, we next decompose the 3-D space into a foliation defined by a spacelike 
vector n. This gives 

(3) i? = le^R^n" + (2 >i? + e n (k AB k AB - k 2 ) 
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Above, e u = — 1, e« = +1. Similarly, one can use a standard expression to write [with q ab = g ab + u a Ub] 

(3) n> (3) n „i/ff 



= Rabcd^e^ (g bd + u b u d ) + K;K up - KK W 
= R ac e^e^ + Raucu&^e-lp) + K^K„p — KK^p 

Using this, we obtain 

^R^n v = n"n u KPK pu - KKp V n»n v + R nn + R nunu 
Putting it all together, we get 

^ ^ R R ~t~ 2 {Ruu Rnn ) ^-Rnunu 

+ f(k) - f(K) - <t>(K) (Bl) 

where /(if) = if 2 — K 2 „, similarly for f(k), and <fi(K) = n^n p (KpK vp — KKp P }. Here, Kp V and /c^s are extrinsic 
curvatures of level surfaces of u embedded in A-D spacctime, and of n embedded in the resultant 3-D space, respec- 
tively This can be rewritten in a better way by defining g^ b = —u a u b + n a n b , and separating out the part depending 
on extrinsic curvatures in II [If, k] = f(k) — f(K) — 4>(K). This yields 

R= ^R + 2R ab g^ b + 2R unun -n[K,k] (B2) 

Note that g ab g^ b = 2, g^ b g^ bc = g^ c , and e ab e cd = 2g^ [d g^ ]b , where e ab = 2u[ a n b ] is the binomial to the surface. 
A slight rearrangement of terms give 

G ab g ±ab = -R ab gi b + R acbd g^ d - <?>R + IL[K,k] 

(B3) 

We stress that the above expression is an identity for any spacelike 2-D surface embedded in a 4-D spacctime, and 
no assumptions such as spherical symmetry or staticity have yet been made. It can be further rewritten as: 

G ab g ±ab = (R uAuB -R nAnB )a AB - ^R + U[K,k] 

(B4) 

where R u a u b = Rabcd.u a u c E b jV) E d B y where E^ are the dyads for the surface. Finally, upon using a AB E^ A ^E h B ^ = 
a ab = g ab + u a u b — n a n b , we arrive at the equation quoted in the text. 



Appendix C: Derivation of Eq. (5) 

Begin with the identity: R^ = — fc°n£ a = V c (fc a V a £ c ) for V[ a fcb] = 0. Now consider the relation fc c V c £° 
-(d z N)k a which is easy to establish for static spacctimcs. We can then prove the following: 

V a (fc c V c n = -(cUV)V -k-k a V a (d z N) 

d 2 N 



-kV • k + 0(z 



2\ 



N 



= - K W-k + 0{z 2 ) - -R± (CI) 

where R± is the Ricciscalar of the t — z part of the metric: — iV 2 dt 2 + dz 2 , and is algebraically equal to —2(d 2 N)/N. 
On the other hand, we also have, from Eq. (3): V ■ k = — In y/a '. 
Putting it all together, we have 

Rtk = K^-lnV^-\R± (C2) 

We now note that, for the 2-D t — n metric, R U nun = — (1/2)jRl [note that u and n form an orthonormaZ basis], 
which finally yields the desired expression 

R$k — Runun = In \fa (C3) 
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Some aspects of this expression have been discussed in the main text. We must also mention that the peculiar 
combination of curvature tensor components above also arises in the analysis in [15] [see Appendices A and B], in 
which area variation of an "acceleration surface" is considered along the observer trajectory [and not along null 
geodesies as in our case] . 
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